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Abstract 

A new family of one-dimensional quantum models is proposed in 
terms of new potentials with a Gaussian asymptotic behavior but ap- 
proaching to the potential of the harmonic oscillator when x — > 0. It 
is shown that, in the energy basis of the harmonic oscillator, the ma- 
trix elements of the Hamiltonian operators of these new models can 
be derived from generating functionals. 
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The progress in algebraic and numerical computation offers one new pos- 
sibilities of analyzing new classical or quantum systems that can not be 
analytically solved. 

The Schrodinger equation with an attractive radial Gaussian potential 
was analyzed with different perturbation methods leading to numerical re- 
sults of convenient accuracy pQ. Other quantum models with Gaussian bar- 
riers were considered in the theory of electronic transport in mesoscopic sys- 
tems [21 [3]. Some time ago we studied a simple one-dimensional relativistic 
model of Gaussian well using generating functionals for finding the matrix ele- 
ments of the Klein-Gordon operator in the energy basis of the non-relativistic 
harmonic oscillator (HO) [I]. In this framework computational methods al- 
lowed us to point out the structure of the discrete energy spectra of the 
Gaussian wells which indicates that inside the well the quantum motion is 
close to an harmonic one. 

Here we would like to continue this study in the case of the non-relativistic 
quantum system, looking for general one- dimensional models with pseudo- 
Gaussian potentials. Our purpose is to define a large family of models whose 
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potentials have a Gaussian asymptotic behavior but approaching to the HO 
potential near x ~ 0. The second objective is to show how the method 
of generating functionals works in this case, helping us to find the matrix 
elements of the Hamiltonian operators of our models in the energy basis of 
HO. 

The one-dimensional HO of mass m and frequency u has the well-known 
Hamiltonian 

h 2 d 2 mtu 2 x 2 
"» = -2^ + — + V " (1) 
where Vb is an arbitrary ground energy. The energy is measured in units of 
e = hw so that we can write 

H = ^N , N = --^ + e + \, (2) 

using the dimensionless coordinate £ = x/a measured in units of a = hj e, 
and denoting A = 2V /e. The number operator N has the eigenvalues 2n + 
A + 1 where n — 0, 1, 2... is the quantum number of the discrete energy levels 
E n = e[n+\{\ + l)}. 

We propose a Gaussian generalization of HO, defining pseudo- Gaussian 
models (PGM) with new number operators 

N = ~ + W(0 , (3) 
in which the dimensionless potentials 

W(t)= (A + E^ 2fc )exp(-/i£ 2 ), (4) 
are supposed to exhibit Taylor expansions, 

w(o = ^+e+o(e +2 ), (5) 

but without terms proportional with £ 4 , £ 6 , £ 2r . It is not difficult to verify 
that this condition is accomplished if we take 

(A/i + k)n k ~ l 
C k - . (6) 

We constructed thus a large family of models, denoted from now by (A,/i) r , 
depending on dimensionless parameters, A G R and /i > 0, and the integer 
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Figure 1: The first six pseudo-Gaussian barriers (A = 8, fi = 0.2) and wells 
(A = -8, fi = 0.2). 

number r = 1,2,... which is called the order of PGM, unlike the genuine 
Gaussian potential that is of the order r = 0. For A > the potentials (jlj) 
are positively defined representing pseudo-Gaussian barriers but when A < 
we can speak about pseudo- Gaussian wells of different profiles (Fig. [Q. 

The PGM we consider here have two principal virtues: their potentials 
(jlj) have a Gaussian asymptotic behavior (provided fi > 0) and, in the same 
time, the form of these potentials in a neighborhood of £ = is similar to 
that of the HO potential. Moreover, when the order r increases, the potential 
(jlj) closes to the HO potential since, according to Eq. (jSJ), near £ ~ these 
potentials differ among themselves only by terms of the order 0(£ 2r+2 ). An 
important consequence of this property is that 

Km = A + £ 2 , V(GK. (7) 

In other respects, we observe that the PGM becomes HO in the limit of 
\i — > too, when we find again that W(£) — > A + £ 2 . 

The Hamiltonian operators of PGM can be written in terms of physical 
quantities using the parameter e as in the case of H . Thus we have 
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where the physical potential, 

n»>=Hf), (9) 

depends on the physical coordinate x measured in units of a defined before. 
The structure of the energy spectra is determined by the concrete form of 
the potential (Q. Since all the models (A,yu) r have potentials that satisfy 

Voo = lim V(x) = (10) 

\x\— >oo 

we can conclude that all of them have continuous energy spectra in the 
domain S C (H) = [0, oo). However, the pseudo-Gaussian wells with A < 
give rise, in addition, to discrete energy spectra in domains Sd(H) = [V m i n , 0) 
with V m in = Vo = ^eA<0. Obviously, the discrete spectra can have at most 
a finite number of energy levels. 

In general, the energy eigenf unctions of PGM as well as the discrete en- 
ergy levels of the pseudo- Gaussian wells can not be calculated using analytical 
methods. Therefore, we must look for appropriate perturbation procedures 
leading to efficient symbolic calculations on computers combined with satis- 
factory convergent numerical methods. 

In the case of our PGM it is convenient to use perturbations in the en- 
ergy basis of HO, {\n) \ n = 0, 1, 2...}, adopting the technique of generating 
functionals related to the generating functions [5] 

F r (0= Q 1/4 exp(-! + 2£T-T 2 ) (n) 

which yield the HO eigenf unctions normalized in the £-scale as, 

Mn (£) = (£| n ) = _^ = ^Ii^l (12 ) 
V ; ^ ' dr n |T =o V ; 

The matrix elements of an operator X can be derived from the corresponding 
generating functional, 

Z ff)T [X] = JdZF a (0[XF T }(0, (13) 

according to the rule 

(m\X\n) = 1 d^dSZrAX] ■ (14) 

V ml nl 2 m+n \a=r=o 
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In general, the matrix (X) with the elements (j!4p is countable since the 
energy basis of HO is countable. The advantage of this method is that for 
HO or PGM the integral ffl3|) reduces to some well-known Gaussian integrals. 
For example, in the simplest case of HO we obtain the functional 



Z<r, T [No] = (1 + A + 4ar) exp(2ar) 

giving rise to the diagonal matrix elements (m\N \n) = (2n + A + l)5 n 
The generating functional of an arbitrary model (A, /x) r , 



(15) 



Z a , T [N] = Z a , T 



+ Z a>r [W(£)] 



(16) 



can be also calculated in terms of Gaussian integrals. Starting with two 
important terms, 



Z a ,T 



Z UjT exp(-//^ 



and the identity ^ 2k exp(— /z£ 2 
result 



" d 2 ' 













exp(2crr) 



exp 



2ar 



(a + r) 



(17) 
(18) 



(— l) k d k exp(— /x£ 2 ), we obtain the final 
Z a , T [N] = Z CT , r [iV]exp(2ar). (19) 



where 

Z a , T [N\ 



- (a-T 



+ 



X + J2(-l) k C k d k 

k=l 



exp 



/i + 1 



(a + T) 



,(20) 



depends on the coefficients (EJ). Now a nice exercise is to show directly that 
lim Z ff , T [N] = KmZ atT [N] = Z^ T [N ] . (21) 

The method of generating functionals offers one the opportunity of using 
computers for finding the matrix elements of the operators N in the energy 
basis of HO. The first step is to derive the generating functional (fl~9l) for 
a PGM of a given order r calculating the function (T201) under Maple or 
Mathematica. Furthermore, one must choose a concrete model (\,fi) r , with 
fixed parameters, for which one has to calculate on computer the matrix 
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elements of N according to the general rule (I14p . Thus one obtains a finite- 
dimensional block (N) representing a truncation of the matrix (N) which 
is countable. Finally, this block can be used for extracting physical results 
(e. g. the energy levels in pseudo- Gaussian wells). We note that in this 
approach the truncation of the matrix (N) takes over the role of the standard 
perturbation procedure. Therefore, increasing the dimension of the block (N) 
one may improve the accuracy of results. 

Our preliminary tests on computers show that the numerical results ob- 
tained for our PGM are satisfactory convergent. Moreover, interesting be- 
haviors can be observed even in the case of the models of lowest orders. 
Thus, it seems that the first discrete energy levels of the pseudo- Gaussian 
wells tend to respect the rules of HO while at the limit of separation between 
the discrete and continuous spectra we meet again the specific effect reported 
in [1], namely an inflexion point of the function giving the dependence of the 
eigenvalues of (N) on the number of level. 

Finally, we conclude that the potentials proposed here realize a gradu- 
ate balance (depending on r) between the short-range HO dynamics and the 
Gaussian asymptotic behavior. We hope that these models should be useful 
for designing new quantum electronic devices in semiconductor heterostruc- 
tures. 
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